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The majority of quantum simulators treat simplified one-band strongly correlated models, whereas 
multiple bands are needed to describe materials with intermediate correlation. We investigate the 
sensitivity of multiband quantum simulators to: (1) the form of optical lattices (2) the interactions 
between electron analogues. Since the kinetic energy terms of electron analogues in a quantum 
simulator and electrons in a solid are identical, by examining both periodic potential and interac¬ 
tion we explore the full problem of many-band quantum simulators within the Born-Oppenheimer 
approximation. Density functional calculations show that bandstructure is highly sensitive to the 
form of optical lattice, and it is necessary to go beyond sinusoidal potentials to ensure that the 
bands closest to the Fermi surface are similar to those in real materials. Analysis of several electron 
analogue types finds that dressed Rydberg atoms (DRAs) have promising interactions for multi 
band quantum simulation. DRA properties can be chosen so that interaction matrices approximate 
those in real systems and decoherence effects are controlled, albeit with parameters at the edge of 
currently available technology. We conclude that multiband quantum simulators implemented using 
the principles established here could provide insight into the complex processes in real materials. 
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I. INTRODUCTION 

Systems of cold atoms have been highly successful as 
quantum simulators for simplified one-band models of 
condensed matter, such as the Hubbard model, where 
atoms move in sinusoidal optical lattices [Ij. Real ma¬ 
terials, on the other hand, consist of multiple electronic 
bands, interacting via Coulomb repulsion and moving in 
a periodic lattice potential with form JT — k e e 2 Z/\r— rq|, 
where Z is the atomic number, r are the locations 
of electrons, r, are the locations of atomic nuclei and 
k e = I/Attcq. It is therefore of interest to examine to 
what extent such a scenario could be reproduced with 
quantum simulators. Painted potentials and holograms 
could be used to make more realistic optical lattices from 
the condensed matter perspective. However, it is not 
possible to reproduce the precise form of the nuclear po¬ 
tential due to optical resolution effects, and this is likely 
to affect the band structure. Another issue is the inter¬ 
action potential, which needs to be of a specific size and 
form. The relative energies of the bands could sensitively 
depend on the form of the external potential and mag¬ 
nitude of interactions between the atoms in the lattice. 
Often states far from the Fermi surface are less impor¬ 
tant, so simulation of solid state phenomena requires as 
a minimum that the bands near the Fermi surface are 
realistically represented. 

Since optical lattices are fixed, they naturally neglect 
phonon effects and are therefore exactly described within 
the Born-Oppenheimer approximation (BOA). The full 
many-body Hamiltonian of a solid state system within 
the BOA contains only three terms: Kinetic energy of 
the electrons, interaction potential between electrons and 
lattice and the interaction between electrons. Since the 
kinetic energy terms of electrons and atoms are the same, 
we can carry out a complete examination of implemen¬ 


tation strategies for quantum simulators of many-body 
condensed matter Hamiltonians within the BOA by in¬ 
vestigating two key aspects: 

1. Establishing to what extent fundamental limita¬ 
tions in the form of optical lattices affects the band 
structure of quantum simulators that emulate the 
multi-band physics of materials. To the best of our 
knowledge this is the first time such an analysis has 
been carried out. We examine the use of painted 
potentials or holograms to form optical lattices that 
are good approximations to the periodic lattice po¬ 
tential. This can be achieved by convolving the 
— k e e 2 Z/r potential (subject to an energy scaling) 
with a Gaussian beam, such that when the beam 
waist is small compared to the optical lattice spac¬ 
ing, the —k e e 2 Z/r potential naturally emerges and 
when beam waist is large the sinusoidal lattice is 
recovered. 

2. Identifying the best electron analogues to use in 
such optical lattice systems. None of the avail¬ 
able electron analogues is ideal. We describe the 
strengths and limitations of systems of cold ions, 
dressed Rydberg atoms, polar molecules and cold 
atoms interacting via Feshbacli resonance for ap¬ 
proximating the interaction between electrons in 
condensed matter systems. We also discuss tem¬ 
perature constraints and decoherence. 

This paper uses graphene and BN as toy systems with 
well understood electronic structures that can be used to 
assess the proposed quantum simulation schemes. How¬ 
ever, the general principles established here can be used 
to guide the implementation of quantum simulators for 
more complex materials. This would provide clean and 
tuneable systems to help understand the physics of con¬ 
densed matter systems in a controlled way. Also, schemes 
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following these principles could be used as much needed 
test beds for condensed matter numerics. 

Before continuing, we briefly describe the various lim¬ 
iting behaviors of condensed matter systems, where the 
ratio of the width of electronic bands to Coulomb inter¬ 
actions gives a rough estimate of the level of quantum 
correlation. Typical cold atom quantum simulators to 
date have concentrated on the strongly correlated sce¬ 
nario of the Hubbard model, which while very interest¬ 
ing from the point of view of experimental development, 
is already well studied using numerical techniques, and 
convergence between these techniques is emerging [2] . In 
this strong correlation limit, i.e. in materials where the 
electrons are well localized to atoms, and electronic bands 
are very narrow, the Hubbard model is accurate. Interac¬ 
tion in this limit is characterized by a single parameter, 
U. The weak U limit is often probed in both optical 
lattice experiments and in numerical simulations. How¬ 
ever, from the point of view of condensed matter systems, 
the single-band Hubbard model is strictly invalid for any 
real system if U is not strong relative to the electron 
hopping between lattice sites, t [3j (perhaps with the ex¬ 
ception of hydrogen under high pressure). Notwithstand¬ 
ing, there is some very interesting physics associated with 
the Hubbard model, including the Mott metal-insulator 
transition [3] (superfluid-insulator transition if bosons are 
used 0), magnetism [B] and potentially superconductiv¬ 
ity (see e.g. [7]). 

The single-band Hubbard model viewpoint is likely to 
miss much of the rich physics found in real condensed 
matter systems. The relative spacing and width of the 
bands is critically important to the size of interband in¬ 
teractions and in many systems the number of interact¬ 
ing bands is large. For systems where the band width is 
small relative to Coulomb repulsion, multiple Hubbard 
parameters can be introduced to maintain accuracy (see 
e.g. [8]). In systems where bands are broad compared 
to interaction energies, correlations are weak and Hub¬ 
bard like approximations break down. The physics in 
the weak correlation limit is quite different and ab-initio 
techniques such as density functional theory (DFT) are 
accurate and popular [5]. 

In the parameter space lying between weakly corre¬ 
lated systems and multi-band Hubbard models, there is 
a region of intermediate correlation where similar energy 
scales are associated with the width of electronic bands 
and the size of the Coulomb repulsion between electrons. 
This regime is difficult to simulate numerically to a good 
level of accuracy, and quantum simulators that can probe 
the physics of such intermediate correlation in systems 
with multiple bands (at least to some degree of approxi¬ 
mation) could offer several benefits: 

1. Cold atom quantum simulators that have multiple 
electronic bands could be used to examine inter¬ 
mediate correlation physics in a controllable way 
that would be useful as a benchmark for numeri¬ 
cal methods, potentially leading to improvements 
in such approaches. 


2. While detailed predictions might be out of the reach 
of quantum simulators, determining the best possi¬ 
ble form and magnitude for the interactions could 
lead to better probes for the rich physics of con¬ 
densed matter systems. 

3. Well separated core states are important for ensur¬ 
ing that the states near to the Fermi surface have 
the correct symmetries, since all states in a 1 band 
optical lattice have an s character, whereas in con¬ 
densed matter systems they may have p, d or / 
symmetry. 

4. Also, interband interactions between itinerant and 
core states can be very important in some cases, 
such as the heavy fermion materials. 

Typical experimental implementations of cold atom 
quantum simulators have focused on the Hubbard model 
FToUn) , and other simulators made from cold ions or Ry¬ 
dberg atoms have been constructed to examine spin sys¬ 
tems IHHIB3- In these systems, atoms or ions in the sim¬ 
ulator represent electrons in a simplified material, and 
an optical lattice represents the external potential from 
a periodic array of nuclei pQ. A wide range of theoreti¬ 
cal studies have focused on the strong correlation limits 
mmm- We also note that quantum simulators of the 
Hubbard and other models can be constructed using cold 
ions, by using the resulting spin systems in combination 
with the Jordan-Wigner transformation mm- 

On the other hand, the weak correlation limit has had 
limited attention [2D]- Of particular interest is the ap¬ 
plication of density functional theory to the band struc¬ 
tures of hydrogen like cold atom systems interacting via 
(5-function like potentials in a shallow sinusoidal opti¬ 
cal lattice, simulated by Ma et al. [21). The all atom 
(all electron) methods used here and in Ref. [21] work 
in continuous space, and should not be confused with 
methods such as Betlie-ansatz local density approxima¬ 
tion (BALDA) which focus on Hubbard models on dis¬ 
crete lattices [22]. Before continuing, we wish to make it 
clear that we are not simulating Hubbard models in this 
paper, rather we are applying density functional theory 
techniques to cold atoms moving in continuous lattices. 

There have been a range of studies where optical lattice 
experiments have moved away from simple square and cu¬ 
bic lattices to investigate more specialized systems. For 
example, honeycomb lattices typical to graphene have 
been implemented (23] [24j. The difference between these 
approaches, and the system proposed here is that the for¬ 
mer simulate a single correlated band, whereas we pro¬ 
pose quantum simulators that can treat all bands of a 
multiband system, including core and itinerant bands, in 
a single simulator. This is far closer to the physics of a 
real solid, and would make it possible to study how tun¬ 
ing the interactions in a material change its properties. 
Other benefits over single band quantum simulators were 
discussed above. We use graphene and BN here because 
they are well understood systems, but the aim is that 
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the ideas presented here can be applied to other, more 
complex, materials. 

This paper is structured as follows. In Sec. [TT]we dis¬ 
cuss the full electron Hamiltonian of a condensed matter 
system, and the considerations for simulating it using an 
optical lattice system. Ideally an optical lattice for con¬ 
densed matter systems would faithfully reproduce the pe¬ 
riodic — k e e 2 Z/r potentials due to the lattice, which im¬ 
poses stringent conditions upon the design of the optical 
system. We also describe the density functional approach 
used to study the formation of energy bands in systems 
with an approximate lattice potential due to an experi¬ 
mentally realizable optical system. In Sec. [Hi] we study 
specific simulated graphene and boron nitride monolayer 
systems. In Sec. [IV]we discuss the magnitude and the 
form of interactions between different electron analogues, 
with the aim of establishing the best system to use in 
experimental implementations. In Sec. [V]we discuss op¬ 
erating temperatures and decoherence effects. Finally we 
summarize and conclude in Sec. IVII 

II. MODEL AND METHOD 

In this section, we discuss our model for optical lat¬ 
tices formed by painted potentials and holograms and the 
method we use to investigate it. We start by discussing 
interactions in a typical condensed matter system. 

The well-known full electron Schrodinger equation for 
a condensed matter system has the form, 

n 2 

~ E V E+E ^ext(T'ia)V’+E ^nt (rij)tp = -E’cm'0 

e i ia i<j 

(!) 

where the wavefunction if) = 'ip(ri , ,r n ) is a many- 
body wavefunction depending on all coordinates of n elec¬ 
trons, displacement between different electrons is r, 3 = 
Tj — rand ri a = r a — r is the displacement between 
electrons and nuclei. In the condensed matter system 
fext (ria) = —k e e 2 Z/ri a , m e is the mass of the electron 
and the interaction potential Vi nt (r , i j) = k e e 2 /rij. The 
direct interaction between nuclei is neglected in Eqn. |T| 
and nuclei are static. Such an approximation is known as 
the Born-Oppenheimer approximation, and is common in 
band structure calculations. In the following, a nn is the 
intersite spacing in the optical lattice, a nn ,CM is the in¬ 
teratomic spacing in the condensed matter system to be 
simulated, M is the mass of the cold atoms in the simula¬ 
tor, and m e the mass of the electron. We use the intersite 
spacing instead of the lattice constant since there are 2 
atoms per unit cell in graphene and BN, and blurring 
effects of the optical system will become relevant on the 
intersite distance before they are relevant on the scale of 
the lattice constant. This will also be the case for more 
complicated unit cells. 

In a quantum simulator, interacting cold atoms, DRAs, 
ions or polar molecules with mass M replace electrons, 
and an optical lattice replaces the periodic potential 


formed by nuclei, but the system naturally has a very 
similar Hamiltonian, 

i ia. 

+ E V ^ T \ r ij)^ AT) = Ekt^ AT) ( 2 ) 

i<j 

where V^xE represents the external potential due to the 
optical lattice and V[ rlt is the interaction between the 
electron analogues. In an optical lattice system, the pe¬ 
riodic potentials representing nuclei are genuinely static 
and the Born-Oppenheimer approximation is exact. This 
corresponds to neglect of phonons. Treating the effects of 
phonons with optical lattices is beyond the scope of this 
paper, and has been discussed elsewhere (See e.g. Ref. 
[23 and ( 25 ]) . For cold atoms, DRAs and cold ions, the in¬ 
teraction potential between the electron analogues may 
be significantly different in magnitude or form from the 
Coulomb interaction between electrons in a crystalline 
material. 

The energy scales in the optical system are necessarily 
lower, because the lattice has a larger size (with changes 
from atomic scales ~ 10~ lo m to optical length scales of 
between ~ 10 -7 and 10 _6 m) and the electron analogues 
(the cold atoms, polar molecules or ions) are heavier, 
which can be understood via a simple scaling argument. 

It can be helpful to rewrite the equation by intro¬ 
ducing the dimensionless variable r' , related to r via 
r = a nni cM^- In this form, 

h 2 

-X-2- V E + ^( fl nn,CM r')lp' = E C m^' 

Zm e a nn,CM 

hence, multiplying by a factor m e a 2 n CM and rearranging, 
the equation reads: 

h 2 

- V'V = m e«nn,CM (E ~ V (a^CM*-')) tf>' (3) 

The goal in the quantum simulator is to find a sys¬ 
tem with an identical second order differential equation. 
The optical lattice has a different length scale, such that 
r' AT a nn = tat- Performing the same analysis, 

- E V 'V (AT) = Ma 2 nn (Eat - V^ AT \a nn r’)) V> ,(Ar) 

so, in order for the prefactor on the right hand side of 
the equations to be the same, 

Mal n E AT = m e al nCM E C M (5) 

therefore the eigenvalues in the cold atom system are 
smaller by a factor of: 

Eat __ m e a nn,CM 

EE “ Mal n 1 j 
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and a similar rescaling is needed for the potentials. Any 
rescaling of the wavefunctions cancels on both sides of 
the Schrddinger equation, and is not important to the 
discussion here. 

An important consideration in the quantum simula¬ 
tor is the use of optical lattices to generate periodic po¬ 
tentials for the cold atoms. The typical sinusoidal form 
of lattice potential is very dissimilar from the attractive 
Coulomb potential generated by nuclei. On the other 
hand, holograms mm and painted potentials [223 offer 
the possibility to make custom optical lattices. 

A starting point for reproducing a Coulomb potential 
with optics is the convolution of a —k e e 2 Z/tqm potential 
with a Gaussian beam with waist, w , which has the form 
exp(—2 r 2 /w 2 ) [3D]. This convoluted potential has the 
form, 


3^ext (x cm) = —fc e e 2 Z'erf(rcMV / 2/w;cM)/ r CMi (7) 

such that the —k e e 2 Z/r potential is recovered as w —> 0. 
This expression is written for condensed matter length 
scales. We require that V /Vat = a ? m CM m e/«nn^, 

luAT = ann^CM/ann.CM an( T r AT = CM/0-nn,CM On 

rescaling for optical lattice length scales. Therefore the 
optical lattice potential should be, 

h ^7 ^ nn ’ CM™e erf(r atv/2/ wat) 

KeC al n M a nn ,cMr A T/«„n W 

, 2 7 «nn,CM»»e erf {rAT/W at) , n . 

k e e Z -—-(9) 

«nn M TAT 

The minimum value of w = A/,/2, where is the wave¬ 
length of the light used to paint the optical lattice, i.e. 
the full beam width is ~ \l (for recent experiments for 
addressing individual lattice sites, see e.g. Ref. [31]). 
Thus, in any optical lattice system the sharp features of 
the —k e e 2 Z/r potential are blurred. 

We also investigate the effects of an additional aug¬ 
mentation where the beam intensity is increased at r = 0 
leading to the optical lattice potential, 


e T W)=- 


(AT) = 2 Z 

ext a nn M 


erf(r-\/2/'ic) 


D 


exp 


2r 2 

w 2 


m 

which gives a slightly better approximation to the 
Coulomb potential if the dimensionless constant, D , is 
not too large. Note that there is an a nn .cm ni e /a n n M 
factor associated with the second term, since w also 
scales. A comparison between these potentials and the 
ideal Coulomb potential can be seen in Fig. [Tk. The 
curve marked D = 0 shows the convolution between 
the Coulomb potential and Gaussian beam. As D is 
increased up to around D = 1.2, the range over which 
the exact Coulomb potential (labeled 1/r) matches with 
Eq. 10 increases significantly. This allows for far bet¬ 
ter agreement between lattice potentials without any re¬ 
quirement for smaller wavelength beams, and is a critical 
observation for the arguments presented in this paper. 


The nuclei in condensed matter systems are arranged in 
periodic arrays, and sums of these potentials can be seen 
in Fig. Eh- These have a Coulomb form for small w, 
and eventually gain a sinusoidal form as w is increased. 
This can be more clearly seen in the Fourier coefficients 
of the array (Fig. EJ-)> which could be used to identify 
simplified optical lattice forms for cases where w is large 
enough that only a small number of spatial harmonics 
are needed to represent the external potential. 

A similar argument requires the ratio of any energy 
scales associated with interactions between the electron 
analogues to be m e a^ n CM /Ma 2 n . We will discuss inter¬ 
actions later in the paper. 

As previously noted, there are two aspects to this pa¬ 
per. The first examines the effects of optical resolution 
on the band structure, making use of the local density ap¬ 
proximation (LDA) of density functional theory (DFT). 
Specifically, we use a modified version of the ELK imple¬ 
mentation [33 ] to examine the band structures and par¬ 
ticle density when the external potential has the form in 
Eqn. [10} ELK is an all-electron linear augmented plane 
wave (LAPW) DFT code. For simplicity and necessity 
in the DFT calculations, we approximate the interaction 
between electron analogues with a k e e 2 /r form. This is 
accurate for cold ions. For cold Rydberg atoms or po¬ 
lar molecules, this gives us a rough method for analyzing 
the effects of the optical lattice resolution, but not the 
effect of non-Coulomb interactions on the band struc¬ 
ture. Choosing a different form of interaction would in¬ 
volve the extremely complicated process of recomputing 
and parameterizing exchange-correlation potentials. The 
second aspect involving the calculation of dipole-dipole 
interaction strengths between dressed Rydberg atoms is 
carried out using van Vleck perturbation theory via a 
custom Mathematica script. In section |IV| we discuss 
the ways in which differences in the interaction matrices 
between Coulomb and dipole-dipole interactions can be 
minimized. DFT calculations are used to examine the ef¬ 
fects of changing interaction strength on the order of the 
difference between the Coulomb and Rydberg interaction 
matrices. 


III. EFFECT OF OPTICAL LATTICE 
RESOLUTION 


Calculations are made for low dimensional quantum 
simulators for materials, since 2D lattices are the sim¬ 
plest to form using painted potentials or holograms. The 
two systems selected for study are a graphene simulator 
and a simulator for the closely related boron nitride, since 
both materials have two dimensional honeycomb lattices 
and weak electronic correlation. Numerical simulations 
of the quantum simulator are made with increasing waist 
sizes w, until the band structure near the Fermi surface 
fundamentally changes form. In the following, results la¬ 
beled w = 0 represent the full electron system calculated 
with the unmodified ELK code (i.e. the full condensed 
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FIG. 1. (Color online) Panel (a) shows an example of a sin¬ 
gle Coulomb potential convoluted with a Gaussian of waist 
w = a nn , and the same potential augmented with an addi¬ 
tional Gaussian beam at r = 0, which makes the potential 
a better approximation to the ideal —Zk e e 2 /r form. Panel 
(b) shows a ID periodic array of potentials, augmented with 
a Gaussian at r = 0 with D = 1.24. Panel (c) shows the 
Fourier cosine coefficients, Fx of a ID array of potentials with 
D = 0 and D = 1.24, to demonstrate how similar features in 
the optical lattices can be established using higher order spa¬ 
tial harmonics, where X is the number of the coefficient. All 
quantities plotted are dimensionless. 


matter systems with an array of —k e e 2 Z/r nuclear po¬ 
tentials). 

Figure [2] shows the predicted bandstructure of the 
graphene simulator for a set of different w. Results us¬ 
ing the full electron local density approximation (LDA) 
of DFT calculated using the unmodified ELK code are 
shown towards the top left labeled as w = 0. As the waist 
is switched on, the band structure relating to the core 
states is modified, and the core states increase in energy. 
As w changes, the band structure maintains a similar 
structure close to the Fermi surface until w ~ 0.149a nn , 


where a band below the Fermi surface exceeds the en¬ 
ergy of the 7r bands at the M point, and heads to¬ 
wards the Fermi energy. Between w = 0.164a nn and 
w = 0.171a nn , this band hits the Fermi surface, with 
the zero band gap semiconductor becoming a semi-metal. 
The band structure near the Fermi surface maintains its 
form up to higher values of w when the Gaussian cor¬ 
rection to the optical lattice potential is introduced. For 
D = 1.24, the physics at the Fermi surface changes be¬ 
tween w = 0.231ct nn and w = 0.239a nn . It may be possi¬ 
ble to use higher spatial harmonics to build up the opti¬ 
cal lattice in place of painted potentials, reference to Fig. 
[ljc) indicates that for w > 0.23a nn the second and third 
spatial harmonics of the Fourier series are important for 
obtaining the correct band structure, but that higher or¬ 
der spatial harmonics contribute only ~ 1.5% of the lat¬ 
tice shape and might be neglected. The relevance of this 
point to operating temperatures will be considered later 
in the paper. Since the core energies become higher than 
valence state energies in the limit that the lattice be¬ 
comes sinusoidal, it is not possible to make a multiband 
quantum simulator simply by adding more atoms per site 
to the modified sinusoidal lattices described in e.g. Refs. 
mi Ml- This demonstrates an advantage to using the 
more realistic optical lattice potentials proposed in this 
paper. 

As a zero bandgap semiconductor, graphene has states 
at the Fermi surface, but it is also instructive to examine 
a band insulator. Figure [3] shows the BN bandstructure 
for several waist sizes (again, the results from the unmod¬ 
ified ELK code are shown, labeled as w = 0). The bands 
start overlapping incorrectly at around w = 0.124a nn - 
The effects are particularly pronounced at the M point, 
and in the case of BN, the band structure goes from hav¬ 
ing a direct gap at the K point, to an indirect gap between 
M and K points. 

Figure [4] shows the electron analogue density for ex¬ 
treme values of w for the graphene quantum simulators. 
The atom (electron) density is much more spread out for 
larger waist sizes. When the additional Gaussian beam at 
r = 0 is added, the particle density becomes more local¬ 
ized, but the band structure at the Fermi surface (Figs. [2] 
and [3]) breaks down when the particle density around the 
nuclei obtain a similar size. This can also be seen for the 
BN analogue in Fig. [5} It is clear that having an appro¬ 
priately localized electron analogue density is key to ob¬ 
taining the correct band structure near the Fermi surface, 
however, it is also important that the lattice potential is 
not significantly deeper than —k e a nn cMm e e 2 Z/a nn Mr 
(which would also localize the electron analogues), as this 
can also lead to significant changes in the band structure. 


IV. INTERACTIONS BETWEEN ELECTRON 
ANALOGUES 


The analysis in the previous section has implications 
for the smallest size of the intersite distance, and due to 








6 



r m r kmk r m r kmk r 
k [ a nn 'l k Kn" J 


FIG. 2. (Color online) Electronic band structure of a quan¬ 
tum simulator for graphene with a variety of feature widths. 
For comparison, results from the unmodified ELK code are 
labeled w = 0 (top left). The Gaussian waist modifies the 
band structure of the core states, and they increase in en¬ 
ergy. Features close to the Fermi surface are well preserved 
until w ~ 0.149o nn , where a band starts increasing in en¬ 
ergy at the M point, and heads towards the Fermi energy. 
Between w = 0.164a nn and w = 0.171a nn , this band hits 
the Fermi surface, with the zero band gap semiconductor be¬ 
coming a semi-metal. When D is turned on a bigger waist 
size can be used. For D = 1.24, the physics at the Fermi 
surface changes for w > 0.239o nn . For the condensed mat¬ 
ter system with w = 0, energy has units of Hartrees (Ha) 
and for the optical lattice systems (w 0) energy has units 
of {m e a^ n CM /Mann) Hartree, which we will denote as scaled 
Hartrees (S. Ha). 1 Hartree « 27.21eV. This convention is 
used throughout, to allow direct comparison between optical 
lattice and condensed matter systems. Ek is quasi-particle 
energy, fc is wavenumber and high symmetry points (F, M 
and K) are defined in the usual way [? ]. 
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FIG. 3. (Color online) Electronic structure of a quantum 
simulator for BN with a variety of feature widths. Again, the 
results from the unmodified ELK code are shown, labeled as 
w = 0. The bands start overlapping incorrectly at around 
ui/a nn = 0.124. The effects are particularly pronounced at 
the M point, and in the case of BN, the band structure goes 
from having a direct gap and the K point, to an indirect gap 
between M and K points. Quantities are defined as in Fig. [2] 


the energy scaling associated with change of length scales 
will affect the size of interaction that is required between 
electron analogues. The smaller the lattice constant, the 
larger the energy scaling factor m e a^ nCM /Ma^ n and thus 
the higher the energy scales (and operating temperature). 
However, the required interaction strength also increases 
as length scales decrease, which could increase decoher¬ 
ence (see Sec. [V]). This indicates a trade-off between 
achievable interaction strengths and operating tempera¬ 
tures. Before we continue, we estimate the shortest pos¬ 
sible distance between lattice sites. 

As we have previously stated, w = Al/ 2. According to 
Sec. Ill the smallest distance between lattice sites that 
can be used to reproduce the physics of the Fermi surface 
is found from w = 0.231a nn ,best for a graphene simulator 
when D = 1.24, corresponding to a nn ,best = A_l/(2 x 
0.231) = 2.16Al. For a blue-green laser, A_l = 532nm 
corresponding to a nnj b es t = 115/um. We will use this 
value of X L throughout. 

Examination of the spatial harmonics of the Fourier 
series indicates that it might be possible to use lattices 
with even smaller constants. While we are unable to 
directly analyze sinusoidal lattices using the ELK code, 
the Fourier analysis in Fig. [TJc) indicates that only the 
fundamental mode plus two spatial harmonics is needed 
to reproduce the optical lattice for w ~ 0.23a nn to around 
1.5% accuracy. 

It would be exceptionally challenging to build the 
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FIG. 4. (Color online) Charge density for graphene quantum 
simulator with (a) w/a nn = 0 (b) w/a n n = 0.149. Panels (c) 
and (d) show the effect of adding Gaussian corrections with 
D = 1.24 to the lattice potential when w/a nn = 0.149 and 
w/a nn = 0.224. Inclusion of the Gaussian correction via D 
localizes atoms to lattice sites, leading to a better approxima¬ 
tion to the condensed matter system. 



FIG. 5. (Color online) Charge density for BN quantum sim¬ 
ulator with (a) w/a nn = 0 (b) w/a nn = 0.109. 


Fourier series using a superposition of standing waves 
formed with lasers with different wavelengths. However, 
it might be possible to use holographic potentials to in¬ 
troduce multiple scales. For the three spatial harmonics 
needed to reproduce the lattice, the fundamental mode 
would be chosen to have lattice constant 3 Al/ 2, the sec¬ 
ond harmonic A l and the third harmonic to have the 
smallest possible spacing of Ar./2. j??i 

We therefore take a lattice constant of d nnj b es t = 3A^/2 
as an absolute limit for the optical lattice built from spa¬ 


tial harmonics (which is slightly smaller than the painted 
potential system). This corresponds to d nni best = 798nm 
using a blue-green laser [Mj. We present results for both 
lattice sizes, and the conclusions of this paper do not rely 
on this smaller lattice constant. Rather, the implemen¬ 
tation of smaller lattices using Fourier components may 
provide an interesting route for improving experimental 
implementations. 

The interaction between the electron analogues must 
be tuneable and close to the values required by the scaling 
argument of Sec. [TT] Here, we mainly focus on dressed 
Rydberg atoms and briefly discuss the relative advan¬ 
tages and disadvantages of systems of polar molecules, 
cold atoms interacting via Feshbach resonances and cold 
ions. For the purpose of calculating the scaling factors, 
the interatomic spacing in graphene is a nnj cM = 1.42A. 

The use of non-Coulomb interactions between electron 
analogues depends critically on the sensitivity of the band 
structure to the interaction strength. We can estimate 
the modification of the band structure in response to 
changes in the interaction potential using DFT, which 
is shown in Fig. [G] A linear dispersion around the K 
points can be found for interactions between V' = 0.8R 
and V' = 1.9R, a scaling factor a little over 2. Here, 
V is the unsealed interaction strength. One of the main 
effects of increased interaction is an increase in the en¬ 
ergy of the core band, and a decrease in the widths of the 
bands, consistent with greater localization of the electron 
analogues. This indicates that the qualitative features of 
the band structure could emerge in a multi-band quan¬ 
tum simulator even if the interaction strengths are only 
approximately similar. 


A. Dressed Rydberg atoms 

In this section, we compare the size of the interaction 
matrices between DRAs with those for the Coulomb in¬ 
teraction. Rydberg atoms interact by long range van 
der Waals forces. By tuning a laser close to the ground 
state to Rydberg state transition, we can make a super¬ 
position of the ground and Rydberg states, where the 
amplitude of the Rydberg state, a, is kept small. In 
this case, the mostly ground state atoms can interact via 
the Rydberg states. If the Rabi frequency of the atom- 
laser interaction is fl 2 ph and the detuning of the laser 
from the ground-Rydberg state transition is A 2 P h, then 
a = fl2ph/2A2 P h- This means the interaction is laser 
controlled, and we benefit by suppressing the effect of 
the finite Rydberg state lifetime, since the lifetime of the 
superposition state is longer by a factor of 1/a 2 , i.e. the 
probability of finding the atom in the Rydberg state is a 
factor a 2 smaller. In order to assess the interaction be¬ 
tween DRAs, we have chosen 43 Ca and 87 Sr which have 
4s 2 and 5s 2 ground states respectively, because they are 
fermionic, have readily available data, have reasonable 
C coefficients (which we compute following Ref. [35]), 
can be excited to 4snR y dS or 5snR y dS Rydberg states 
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FIG. 6. Effect of scaling the interaction on the band structure 
of simulated graphene. Here, w = 0.0745a nn . The K points 
are well defined between V' = 0.8V and V' = 1.9V, a scaling 
factor a little over 2, with the separate Is core state separated 
from the n and sp 2 bands until V' = 1.5V. The main effect of 
increased interaction is an increase in the energy of the core 
band, and a decrease in bandwidth size indicating greater 
localization of the free particles. Quantities are defined as in 

Fig- 0 


with nearly spherically symmetric interactions and be¬ 
cause there is exists a transition with extremely long life¬ 
time which greatly suppresses heating in the Rydberg 
dressing scheme (see Sec. [vf. We note that lighter al¬ 
kaline earth metals are likely to have similar transitions 
and narrow linewidths and might have other advantages 
(e.g. higher operating temperature). 

For almost spherically symmetric interactions we wish 
to excite the 4sn Ryd s or 5sn Ryd s Rydberg states respec¬ 
tively, where ?r Ryd is the principal quantum number of the 
Rydberg state, and this means the laser dressing must be 
via two photons. The first photon is from a laser tuned 
close to the narrow linewidth 1 S’ 0 —>- 3 Pi “intercombina¬ 
tion line” with Rabi frequency uj\ and detuning <5i (for 


43 Ca the linewidth is 27 t x 6.5 kHz and wavelength 657 
nm; for 87 Sr the linewidth is 27r x 7.5 kHz and wave¬ 
length 689 nm). The second photon is from a laser tuned 
close to the 43 Ca 4s4p —>• 4sn Ryd s transition (or the 87 Sr 
5s5p —> 5sn Ryd s transition), with Rabi frequency oj 2 and 
detuning <5 2 . The effective “two photon” Rabi frequency 
is then fl 2 ph = Wiu; 2 /2(y 1 , and the two-photon detuning 
A 2p h = <5i + 62 [31]. A schematic of the dressing scheme 
for 43 Ca is shown in Fig. [ 7 ] The 87 Sr dressing scheme 
is similar. The interactions between atoms are fully de¬ 
termined from A 2p h and fl 2p h, however the values of 5i 
and iOi can affect heating (see Sec. |V|. The interactions 
between dressed Rydberg atoms are highly tuneable and 
with the right dressing scheme, almost spherically sym¬ 
metric interactions between the dipole moments of the 
atoms can be induced m- 

For simplicity, we consider the interactions of two 
atoms in the presence of the dressing light. The dressed 
Rydberg atom interaction potential can be calculated us¬ 
ing e.g. van Vleck perturbation theory [35], which gives, 

i/Tot _ ( e± “ 2A 2p h)A 2p hH2 ph + H 4 ph 

± ~ A 2 p h (4A 2 ph -Q 2 ph -2 e± A 2ph )- 1 J 

where e± represent the bare Rydberg-Rydberg interac¬ 
tion (see below). In the limit where e± —>• 0, V± ot —► 
Vl = (^ 2p h — 2A2 ph fl2 P h)/(4A 2ph — A 2p h^| p h)) which is 
due only to the laser dressing, but not the interactions 
of the atoms. We subtract Vl from V± ot to obtain the 
effective interaction potential: 


v± = _ e± ^2 P h _ 

(^2 P h — 4A2 ph ^ ^2e±A 2p h + H| ph — 4A2 ph ^ 

( 12 ) 

When the separation between the two atoms is small, 
the shift of the Rydberg levels by the interaction inhibits 
the excitation of the Rydberg levels by the laser, and so 
the interaction saturates at — fl|p h /(8A 3 ph — 2A 2p h^ 2ph ), 
which is ~ H 2p hQ: 3 when A 2p h 3> H 2p h. To calculate e±, 
we must diagonalize the Hamiltonian for the two-atonr 
states, which takes the form of a 2 x 2 matrix in the basis 
|n Ryd s, n Ryd s), |n' Ryd p, n£ yd p): 


( 0 C/r 3 \ 

^ C/r 3 5 Foi ) 


(13) 


where <5 F or = ( E n' RydP + E n£ yA p) ~ 2 E n RydS is the Forster 
defect, and E nRyd i is the energy of n Ryd / state of the 
atom. In the dressed Rydberg scheme, the laser driven 
transition is always much closer to the state with the 
eigen-value e_|_ > 0, which ensures the interactions are 
repulsive. Therefore, in the following, we will always be 
discussing V + when we refer to the interaction potential. 
The C coefficient is determined by the Rydberg states 
involved in the scattering, following Ref. [35|. We need 
to choose an interaction that will mimic the scaled in¬ 
teraction matrix strength, which depends on the optical 
lattice constant and the mass of the atom used in the 





















9 


(a) 


(b) 


matrices for such an interaction, which have the form, 


Si + ^2 


2 


0) 2 


i 

1 

i 
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I1IZ 5 1 + S 2 = A 2|* 
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FIG. 7. Scheme used to dress the atoms with lasers, (a) 
shows the relevant transitions for driving a single 43 Ca or 
8 'Sr atom via a two-photon transition, with Rabi frequen¬ 
cies LOi and detunings Si, where i=l,2 labels the transitions. 

(b) the two photon transition is equivalent to a single tran¬ 
sition driven with Rabi frequency fl 2 P h and detuning A 2 P h- 

(c) For two atoms, the excitation to the 2 atom states has an 
energy spectrum modified by the Rydberg-Rydberg interac¬ 
tion. The two-photon transition dresses the higher energy 2 
Rydberg-atom state, being far detuned from the lower energy 
2 Rydberg atom state. 


simulator. At the short distances over which the interac¬ 
tion potential is significant (up to about 2 microns), the 
potential can be simplified because e+ « C/r 3 when r is 
small. It follows that 

+ (a 2 — l) 2 r 3 + Ca/fl 2p h(a 2 — 1) 1 1 

so V + has the form V + = Vo/(r 3 + r 3 ), with r c = 

(Ca/f^ph^ 2 — l)) 1 / 3 - N.B. Both a 2 — 1 and a are neg¬ 
ative in the following, so r c is positive. 

We proceed to estimate the strengths of the interaction 


U = J d 3 rd 3 r'V(r — r')N(r)N'(r') (15) 

We assume that the electron (electron analogue) density 
is step like, N(r) = N (r < b) and N(r) = 0 otherwise, 
approximating the form of a Wannier function. There¬ 
fore, N = 3/47 rb 3 . For in-band interactions (which are 
largest), N = N'. Interband interactions are smaller. 
Matching of the interactions on the scale of lattice spac¬ 
ing requires that U/U cou i = m e a 2 nCM /Ma 2 n . The ex¬ 
perimental aim would be to get a match within a range 
of a few 10s of percent, so that changes in interaction do 
not strongly change the band structure (see also Fig. |6|. 

By modifying the principal quantum number of the 
Rydberg state, it is possible to modify the range of 
distances over which the interaction matrices for both 
Coulomb and dipole-dipole interactions have roughly the 
same magnitude. We change the strength of the inter¬ 
action by varying A 2 P h, 7iR y d and f^ph- A summary of 
the values used can be found at the top of Table [l| The 
form of the interactions for the different parameter sets 
is shown in Fig. [8j Numerical evaluations of the interac¬ 
tion matrix elements using Mathematica’s Monte Carlo 
integration routines for the four sets of parameters can be 
seen in Fig. [9] The figures also show the rough bounds 
on the interaction matrices determined by the DFT cal¬ 
culations in Fig. [6] 

It can be seen that while the interactions do not have a 
Coulomb form, the interaction matrices have magnitudes 
within the bounds indicated by the DFT calculations over 
a range of Wannier-like function sizes, from radii slightly 
larger than those shown in the figure, down to around 
25% of the largest sizes shown. The reason for the region 
of approximate agreement is that the plateau in the inter¬ 
action at small r means that Coulomb and dipole-dipole 
interactions have a region where they are tangential. An 
experimental protocol would aim to match interaction 
matrices for Wannier-like functions on the size of an inter¬ 
atomic spacing (i.e. those most associated with the Fermi 
surface) by tuning the dipole-dipole interaction through 
the Rabi frequency, detuning parameter, and principal 
quantum number of the Rydberg state to maximize the 
length scales over which the interaction matrices are of 
similar magnitude. This is likely to be problem specific. 

We can speculate about how this might affect differ¬ 
ent states by using hydrogenic wavefunctions as an esti¬ 
mate. The radius of such wavefunctions scales roughly as 
n Wann- Hence riwann = 1 states are approximately 1/4 of 
the size of nwann = 2 states, where nwann is the princi¬ 
pal quantum number representative of the Wannier state 
(not to be confused with the principal quantum num¬ 
ber of the Rydberg state). The results indicate that the 
interaction matrices of materials containing elements in 
the first 2 rows of the periodic table could be approx¬ 
imately represented using the types of DRAs we have 
discussed here, although there may be difficulties with 
quantum simulation of materials containing elements in 
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TABLE I. Parameters for proposed experimental DRA quantum simulator. The upper part of the table summarize the 
parameters of the quantum simulator for the specified atom, and the properties of the transitions that should be driven. The 
last 2 column are for 25 Mg, for which we do not have parameters for the Rydberg levels, so we estimate using C and feor for 
Ca (see text). We select the same Rabi frequencies for both transitions required for the two-photon scheme, so that uj 2 = wi, 
and the detuning for the upper transition satisfies 82 = A 2 p n — 81 . The lower rows of the table relate to decoherence effects and 
are discussed in Sec. [v] We also show the laser power needed to drive the lower 1 S'o —> 3 Pi “intercombination” line transition, 
assuming that the laser is focused into either a 50 x 50 /im area (for 43 Ca and 87 Sr) or 25 x 25 /im for ( 25 Mg). The lowest 
row gives an estimate of the number of hops per simulation, which is calculated from either the DRA lifetime or the heating 
(whichever is smaller). 


Species 


43 Ca 


87 Sr 


25 Mg 

Transition 

4s 2 —l 4s4p —► 4 sn Ryd s 

5s 2 —> 5s5p —> 5sn Ryd s 

3s 2 —¥ 3s3p —> 3sn Ryd s 

U Ryd 

32 

38 

32 

38 

32 

38 

C (MHz /nn 3 ) 

24.7 

50.1 

29.3 

63.3 

Estimate from Ca 

<W (GHz) 

-16.700 

-8.880 

-10.900 

-5.440 

Estimate from Ca 

A 3Ri (nm) 

657 

657 

689 

689 

457 

457 

r 3Pl (Hz) 

6000 

6000 

7460 

7460 

31 

31 

4at3p, MW/cm 2 

2.770 

2.770 

2.980 

2.980 

0.042 

0.042 

Inter-site distance a nn (nm) 

798 

1150 

798 

1150 

798 

1150 

0 2ph (MHz) 

6.17 

2.55 

6.51 

2.55 

7.07 

2.92 

A 2p h (MHz) 

-34.29 

-14.17 

-46.48 

-17.96 

-35.35 

-14.60 

a 

-0.09 

-0.09 

-0.07 

-0.071 

-0.1 

-0.1 

1 /a 2 

123.5 

123.5 

204.1 

198.4 

100 

100 

Effective Lifetime Ryd. aA 2 T Ry( j (ms) 

1.57 

2.63 

2.59 

4.22 

1.27 

2.13 

detuning ns 2 —¥ nsnp: <5i (GHz) 

9.19 

9.19 

9.19 

9.19 

0.50 

0.50 

Rabi freq. ns 2 —> nsnp\ uj 1 (MHz) 

337 

217 

346 

217 

84 

54 

Laser power (mW) 

436 

180 

320 

126 

3900 

1610 

DRA heating rate E Ryd (nK/ms) 

13.0 

5.4 

7.6 

3.0 

5.0 

2.1 

A eg (nm) 

423 

423 

460 

460 

285 

285 

T eg (GHz) 

0.2 

0.2 

0.217 

0.217 

0.5 

0.5 

A l (nm) 

532 

532 

532 

532 

532 

532 

Lat. heating rate, El (nK/ms) 

0.56 

0.39 

0.21 

0.14 

1.13 

0.78 

Total heating rate, E (nK/ms) 

13.56 

5.79 

7.81 

3.14 

6.13 

2.88 

Hopping time r (ms) 

0.20 

0.4 

0.41 

0.85 

0.12 

0.24 

T\o% (nK) 

8.43 

4.05 

4.17 

2.00 

14.4 

7.2 

hops/sim. ss T 10 % / Et or r Ryd /ra 2 

3.11 

1.76 

1.29 

0.75 

10.58 

8.875 


the 3rd and higher rows. This difficulty is partly related 
to heating, which we discuss later, and might be avoided 
by using different Rydberg dressing schemes to the ones 
discussed here. The interactions associated with electron 
analogues in the smaller core states would be underesti¬ 
mated, and probably the energies of these states would 
be lowered. The Rydberg dressing scheme would need 
careful tuning for each problem to be investigated, so re¬ 
sults from such a quantum simulator may be able to give 
useful insight into specific physical processes in materials, 
but would not be predictive. 

B. Polar molecules 

For completeness, we mention an intriguing possibility 
to use polar molecules. The recent formation of polar 
molecules in cold atom systems means that it is possi¬ 


ble to generate customizable interactions with a simi¬ 
lar form to those of the dressed Rydberg states mm- 
Under certain circumstances, they might offer better in¬ 
teraction properties than cold Rydberg atoms. For ex¬ 
ample the form of the interaction can be controlled at 
very short distances, such that it may be possible to ad¬ 
just the interaction with multiple avoided crossings so 
that the regions where the 1 /r and dipole-dipole interac¬ 
tion curves match are longer [39j. Particular benefits of 
polar molecules are narrow linewidths of the transitions 
that lead to long lived excited states and low scattering 
rates (which would lead to low heating). Dipoles and 
therefore interaction strengths are induced by external 
electric fields, which makes them highly controllable, al¬ 
though for very strong fields there could be issues with 
anisotropy. Since the form of the interactions is very sim¬ 
ilar to that for dressed Rydberg atoms, the arguments of 
the previous section regarding the interaction matrices 
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Atomic separation (fim) 



Atomic separation (jura) 


FIG. 8. Interactions between dressed Rydberg states of pairs 
of 43 Ca atoms, and pairs of 87 Sr atoms. Coefficients are as in 
Tabled 

should apply. We do not have access to sufficient infor¬ 
mation to make quantitative predictions regarding the 
use of polar molecules at the current time. 

C. Feshbach resonances and cold atoms 

Interactions between cold atoms can be magnetically 
tuned using Feshbach resonances |T| , such that they have 
a delta function form, Anh 2 a(B)S(r) / M, where 

cl(B) = a bg ^1 - b Af Bq ) (!6) 

where is the resonance width, a bg is the background 
scattering length, and B the magnetic field mm- 
The ratio of the interaction matrices diverges for very 
small Wannier-like functions. However, we suggest that 
use of weak interactions induced by Feshbach resonances 
might be useful to correct the underestimation in the 
ratio of the interaction matrices on short distances when 
using dressed Rydberg states [42] . For large B fields, 
anisotropy in Rydberg states could become a problem if 
making this correction. 




Wannier wavefunction range, b (pm) 


FIG. 9. Ratio of the interaction matrix elements of Rydberg 
atoms to Coulomb interactions calculated using Monte Carlo 
integration for (a) 43 Ca and (b) 87 Sr. The horizontal lines 
show the maximum and minimum bounds on the interaction 
strength determined using DFT. 


D. Cold ions 


The scaling that is required for the interactions be¬ 
tween particles means that although 1/r interaction is 
available for cold ions loaded into an optical lattice, the 
energy scales are wrong for quantum simulation of ma¬ 
terials by several orders of magnitude. The interactions 
must be scaled by a factor m e a nni cM/Ma nn (as for the 
—k e e 2 Z/r nuclear potential, there is a cancellation of 
factors of a nn and a nn ,cM that emerge when scaling r). 
However, the interactions between ions are fixed by their 
charge, so cold ions in optical lattices will have poorly 
representative band structures. This is confirmed by Fig. 
[6j where increasing the size of the potential by only a fac¬ 
tor of 2 dramatically changes the band structure close to 
the Fermi surface. 
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V. OPERATING TEMPERATURE AND 
DECOHERENCE 


For the final section of this article, we will focus on 
decoherence effects in dressed Rydberg atoms. These ef¬ 
fects can be divided into decoherence by atom heating 
due to the optical lattice and laser dressing, and deco¬ 
herence by atom losses. 


A. Heating due to the optical lattice 


A source of decoherence is the absorption of photons 
used to form the optical lattice and to dress Rydberg 
atoms. This problem has been treated in depth in j43j . 
The optical lattice is painted or generated by holograms 
using a laser of wavelength \ L (detuning A) and Rabi 
frequency fi, distinct from those involved in the dress¬ 
ing of the Rydberg states. Energy is absorbed by each 
atom at a rate El = 7o = T eg (fi 2 /4A 2 ) 2 Er, where 
Er = (/i/Al) 2 /2M is the recoil energy, T eg is the natu¬ 
ral linewidth and A eg is the wavelength of the 1 S'o —> 4 Pi 


transition [J3|- In the expression for El , the factor of 
two appearing in front of the recoil energy given above 
corresponds to the average kinetic energy of an atom re¬ 
coiling from absorbing a laser photon and emitting spon¬ 
taneously into the full 47r steradians solid angle [33] , 
The maximum lattice depth can be calculated by tak¬ 
ing the leading term of the Maclaurin expansion of Eqn. 


10} which can be computed using erf(z) ~ 2z/y / 7r as 
tw T) (°) = -a nn ,CMm e fc e e 2 Z(2 3 / 2 /- v /7f + D)/wa nn M. 
This should be matched to the optical lattice parameters 
via Idattice = ft.0 2 /4A which corresponds to the maxi¬ 
mum laser intensity. A = c(l/A p — 1/A eg ), so the lattice 
depth fixes H. We select A l = 532nm for all atomic 
species. This estimated rate of heating due to the lattice 
with D = 1.24 is shown in Tab. |U 


B. Heating due to laser dressing 


Heating associated with the laser dressing is more 
rapid because the scheme only works when the lasers are 
tuned fairly close to the atomic transitions. The first 
requirement that we must satisfy is that the dressing is 
strong enough. The strength of the Rydberg-Rydberg 
dressed intera ction saturates at ~ Ha 3 . As discussed 
in Sec. IV A since we are driving a 4s 2 to 4snR yc jS in 


43 Ca or 5s 2 to 5snR, yc js transition in 87 Sr, we need a two- 
photon laser excitation scheme, which means excitation 
via an intermediate state. The two lasers have Rabi- 
frequencies of u>i and detuning Si, i=l,2. We want to 
achieve a two-photon Rabi frequency of l^ph = ui^/ZSi 
and a two-photon detuning of A 2 P h = <5i + Si- In this 
scheme, heating will arise from driving the lower tran¬ 
sition, whereas driving the upper transition will lead to 
losses from the Rydberg state. This occurs because the 


Rydberg states decay via a cascade, whereas the 4s4p 
state decays back to the 4s 2 ground state. The heating 
rate derived using similar arguments to Refs. @3} H2] 
is I?R y d = (/i 2 /A 2 p i M)(r3p 1 w 2 /4<5 2 ), the first factor be¬ 
ing the kinetic energy of the recoiling atom, the second 
factor is the scattering rate. If we pick io\ = u> 2 , then 
pRyd = (^ 2 /A3p i M)(r3p 1 H 2p h/2(5i). It is because of the 
heating due to the dressing scheme that we have chosen to 
study the Alkaline-Earth elements, rather than the Alkali 
metals common in cold atom experiments. In the case of 
the Alkali atoms, the intermediate states would be the 
nP states closest to the ground state, all of which have 
T ~ 2-7T x 6 MHz linewidths and correspondingly high 
scattering rates. The Alkaline-earth metals Ca and Sr, 
however, benefit from a low-lying state with linewidth, 
r3 Pl , on the order of 27 t x 7kHz, which is narrow enough 
to allow laser cooling to very close to quantum degener¬ 
acy. 


Heating rates due to dressing are summarized in Table 
[jJ p or 43 Ca, calculations assume driving a two photon 
transition exploiting the 657 nrn 1 So —> 3 Pi intercom¬ 
bination transition (T3 p 1 /2tt = 6 kHz). For 87 Sr, this 
transition has X^p 1 = 689 nm and linewidth T3 Pl /2i: = 
7.5 kHz. To suppress the heating rate requires detuning 
the lasers very far from the one-photon resonances, with 
a corresponding increase in Rabi frequency. The smaller 
values of wi and <5i that we have chosen in this example 
are challenging, requiring large laser intensity (around 
130 - 200 mW when focused into 50 x 50 ^m 2 ). We com¬ 
pute the laser power required by noting that in terms of 
the laser intensity Jo, uq = T^/lo/2I sat , where the satu¬ 
ration intensity is I sat = 2ir 2 hr3 Pl l3Al Pi [3S]. The table 
shows that heating rates are on the order of a few nK / 
ms, slightly larger than the heating from the lattice. 


We estimate that it should be possible to reduce heat¬ 
ing further using Mg. We do not have access to the re¬ 
liable data regarding quantum defects for Mg Rydberg 
levels, but we make the observation that given the simi¬ 
larity between the Ca and Sr Rydberg levels, a suitable 
Rydberg level, giving useful interaction parameters, will 
almost certainly exist for Mg. The Mg intercombination 
line is even more highly suppressed than those of Ca and 
Sr, because Mg does not have the strong spin-orbit cou¬ 
pling of the heavier elements. Hence, the Mg 1 S'i —» 3 Pi 
transition linewidth is two orders of magnitude lower, 
r 3Pl = 2-7T x 31 Hz; for this transition, A 3 Pl = 457 nm, 
so while the recoil energy is 3.55 times larger than for 
43 Ca, the heating rate is significantly reduced. With such 
a small linewidth, and correspondingly small saturation 
intensity, the laser power required to drive the intercom¬ 
bination transition is significantly increased, but with the 
advantage of a slower heating rate: see Table [T] We cal¬ 
culate that about 1.6 W of laser power focused into an 
area of 25 x 25 /xm 2 is sufficient to drive the transition. 
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C. Atom losses 

Next, we discuss atom losses. The main loss will be 
due to the finite Rydberg atom lifetime ~ 21 jx sec, which, 
since the probability of exciting the Rydberg state is a 2 , 
can be expected to be enhanced by laser dressing, by a 
factor of a -2 , (recall that the amplitude of the Rydberg 
state is a = fl 2 ph/2A2 P h)-) The lifetime of the Sr 5s20s 
state is about 3.1 fjs @B], and using the n^ yd scaling, 
we estimate that the Sr 5s38s state lifetime will be 21 
/is (we expect this lifetime will be similar in other alka¬ 
line earths, which are estimated in the same way). For 
our parameters, or 2 > 100. This limits the experiments 
to around 1 to 4 ms. Although this is very short, the 
lasers can be tuned to improve this situation. Finally, 
we note that the dressed atom schemes are subject also 
to collective decoherence effects, which are the subject of 
experimental and theoretical work So far, the 

results of experiment of Ref. (48| achieved improved Ry¬ 
dberg lifetimes, however the increase was far below the 
theoretical maximum. 


D. Consequences of decoherence effects 

For graphene, t ~ 3eV, so the hopping energy 
in the quantum simulator is approximately f at = 
?u e a 2 n CM 3eV//iMa 2 n . Typical timescales to hop be¬ 
tween sites are h/2tAT on a tight binding chain (see 
e.g. [50]) which corresponds well to the group velocity 
at the Fermi surface of a half filled chain. Using the 
group velocity at the Fermi surface of graphene, which is 
v g = 3ta/2h, we find that the typical timescales, r, to hop 
between carbon atoms are slightly slower r = 2?i/3tAT- 

The bandwidth of the tt bands in graphene is 6 1 ~ 
18 eV. Rating the quality of the quantum simulator 
depends on how well resolved the bandstructure needs 
to be. To resolve bandstructure features to 10%, cor¬ 
responds to an energy scale of ro e a 2 n CM /Ma 2 n x 1.8eV 
in the cold atom system. The temperature of the cold 
atom system needs to be low enough to resolve these de¬ 
tails, and can be estimated as Tio% = m e a 2 n CM /Ma 2 n x 
1.8eV /kn, where ks is Boltzman’s constant. Finally, 
the temperatures must remain stable for long enough for 
hopping between lattice sites. Assuming that the initial 
temperature of the simulator is much lower than T 10 %, 
then a rough estimate of the number of hops until the 
simulator is too warm to resolve bandstructure features 
to 10% is T 10 %/tE, where E is the total heating rate. 
This estimate will be better for simulators with higher 
operating temperatures (and larger numbers of hops per 
experiment). The number of hops may also be limited by 
Rydberg lifetimes as r^ yd /a 2 T. The approximate num¬ 
ber of hops per experiment taking the smallest of the two 
values is shown in Table |T| This indicates that so long as 
linewidths remain small, lighter atoms and shorter lat¬ 
tices perform better (have more hops) before heating be¬ 
comes an issue, and that detuning should be maximized. 


Many condensed matter processes require only hops to 
and from a neighboring site, so a few hops should be suf¬ 
ficient to bring the simulator into equilibrium in many 
cases, especially if the interaction is turned on adiabati- 
cally. As an extreme example, just 3 hops are required to 
restore equilibrium to a Bose-Hubbard model undergoing 
a quench from a superfluid to Mott insulator state m- 
Our estimates indicate that Mg DRAs should be used. 

VI. CONCLUSIONS 

We have examined some of the challenges involved in 
making a quantum simulator that can emulate multi¬ 
band materials. Optical lattices are fixed and there¬ 
fore naturally follow the Born-Oppenheimer approxi¬ 
mation where phonons are neglected. Within Born- 
Oppenheimer, the full many-body Hamiltonian contains 
three terms: Kinetic energy of the electrons, interaction 
between electrons and the lattice and the interaction be¬ 
tween electrons. Since the kinetic energies of electrons 
in a condensed matter system and cold atoms in a quan¬ 
tum simulator are the same, we focused on ways in which 
the two interaction potentials could be approximated in 
a multiband quantum simulator. 

1. Our first aim was understanding how optical lattice 
resolution changes the bandstructure of a multi¬ 
band quantum simulator. We have performed DFT 
calculations to examine these effects. This is a 
step towards determining the forms of optical lat¬ 
tice suitable for quantum simulators with the ap¬ 
proximate band structures of real materials. Our 
primary conclusions are that: 

• Resolution effects do not modify the band 
structure around the Fermi surface over a wide 
range of beam waist sizes. 

• Core state energies may overlap valence states 
as vu/a nn becomes large and the lattice poten¬ 
tial becomes sinusoidal. Since core and va¬ 
lence states are typically well separated in real 
materials, standard sinusoidal optical lattices 
may not be suitable for many band applica¬ 
tions. 

• Optical resolution effects are minimized by 
making a nn /w large. However a nn cannot be 
increased indefinitely since operating temper¬ 
atures decrease with a nn . This can be par¬ 
tially mitigated using short wavelength lasers 
and lighter atoms to increase operating tem¬ 
peratures. 

2. Our second aim was to identify electron analogues 
with a suitable magnitude of interaction matrices. 
Our second set of conclusions are as follows: 

• DFT calculations indicate that bandstruc- 
tures maintain qualitatively similar forms 
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even if interaction strengths vary by around 
a factor of 2. 

• The strength of interactions between cold ions 
is fixed by the electron charge and is too strong 
for probing multiband physics. 

• Dressed Rydberg atoms have highly tuneable 
dipole-dipole interactions. Calculations using 
van Vleck perturbation theory show that in¬ 
teraction matrices are similar to those of a 
Coulomb potential over a range of Wannier 
function sizes. Quantum simulation is possi¬ 
ble, but limited to a few hops. 

• Polar molecules are also highly tuneable. We 
have not made detailed calculations here, but 
they could offer a similar benefits to Rydberg 
atoms. 

• Cold atoms interacting via Feshbach reso¬ 
nances have interaction matrix elements that 
diverge relative to those of the Coulomb in¬ 
teraction for Wannier functions with small ra¬ 


dius. 

Our interpretation of these results is that simulation of 
multiband materials is possible using DRAs, but it will 
be approximate rather than predictive, and right at the 
limit of currently available Rydberg dressing and optical 
lattice technology. This paper discussed the toy systems 
of graphene and BN because of the need to make con¬ 
tact with the well understood electronic structure of these 
materials. The goal is that more complex systems could 
be simulated. Interactions would need tuning for each 
system to be investigated, and could not be selected ab- 
initio. Schemes following these principles could be used 
as test beds for numerical condensed matter techniques. 
Future work should focus on improving the interactions 
between the electron analogues while decreasing decoher¬ 
ence effects. For example, single phonon Rydberg dress¬ 
ing schemes with very low heating rates may be possible 
using atoms with a P ground state, such as In (523. It 
would also be very interesting to find ways to introduce 
tuneable phonon effects into multiband quantum simula¬ 
tors. 
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